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We review the recent theoretical progress in heavy quarkonium spectroscopy within the boundstate theory based 
on perturbative QCD. New microscopic pictures of the heavy quarkonium systems are obtained. 



1. Introduction 

Recent advent of the boundstate theory based 
on perturbative QCD enabled accurate descrip- 
tions of the nature of the heavy quarkonium sys- 
tems such as bottomonium states. Through the 
progress became more and more clear the impor- 
tance to eliminate properly contributions of in- 
frared (IR) degrees of freedom in the theoretical 
descriptions of these systems. Indeed these sys- 
tems comprise natural IR cutoffs both of spacial 
and temporal dimensions. As a result, once we 
decouple the IR contributions properly, we find 
much better convergence of perturbative expan- 
sions and a good control over the theoretical pre- 
dictions. In order to realize IR decoupling in a 
systematic way, the language of renormalons and 
their cancellations have played central roles: un- 
certainties of the perturbative expansions are es- 
timated from their asymptotic behaviors based on 
renormalon dominance picture. 

In the heavy quarkonium systems, we may clas- 
sify the mechanisms of IR decoupling into two 
categories: (I) The spacial size of a quarkonium 
state acts as an IR cutoff. This leads to cancel- 
lation of 0(A) renormalons, with residual 0(A 3 ) 
renormalons. (II) Offshellness of the heavy quark 
and antiquark acts as an IR cutoff in temporal 
dimension. This leads to cancellation of C(A 3 ) 
renormalons, with residual 0(A 4 ) renormalons. 
These aspects point to new and detailed physical 
pictures of the heavy quarkonium states. 

2. Cancellation of 0(A) Renormalons 

Cancellation of 0(A) renormalons [1] is a re- 
alization of the fact that the system under con- 
sideration is color-singlet and has a spacial size 
R much smaller than the typical hadron size 
Ai?< 1. 



10 


.5 


> 


10 






o 




9 


.5 


>> 




'-.i. 




u 
o 

a 


9 


a 




8 


.5 



Power— law potential 
Cornell potential 

Log potential 




|4 

J/if> V>(2S) 

Fig.l: -Etot(f) and typical phenomenological po- 
tentials. Constants are added to make all curves 
coincide at r = 1 GeV^ 1 . Arrows at the bottom 
show the r.m.s. radii of the heavy quarkonium 
states. The figure is taken from [2]. 

2.1. Total energy of a static bb pair 

Let us first neglect the kinetic energy of b and b 
and examine the total energy of a static bb pair. It 
is defined as the sum of the static QCD potential 
and the pole masses of b and b: 



E tot (r) = 2m biPO \ c + VqcdM- 



(1) 



*Talk given at "International Conference in Quantum 
Chromodynamics (QCD02)", Montpellier, France, July 2— 
9, 2002. 

tTU-667, September 2002. 



The 0(A) renormalons contained in the pole 
mass and the QCD potential are cancelled if we 
rewrite the pole mass in terms of the MS mass, 
nib = m^ ls (m^ IS ). We thus substitute mfc iPO i = 
mb(l + c\ots + C2a 2 s +C3a|) to eq.(l) and expand 
E io t{r; m,b,as(p), in as(fJ-) upto0(a|). Then 
the perturbative expansion becomes much more 
convergent as well as much less dependent on the 
renormalization scale [i (see eqs.(3)-(6) below). 
The remaining renormalon is of order A x (A • r) 2 . 
Fig.l shows E tot (r) corresponding to the present 
values of the strong coupling constant. They 
agree well with typical potentials used in phe- 
nomenological model analyses within the uncer- 
tainty estimated from the remaining 

A 3 r 2 

renor- 
malon (indicated by error bars), in the range rele- 
vant to the bottomonium and charmonium states. 
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Recently several comparisons have been made 
[3,4,2,5] among the perturbative predictions of 
the QCD potential (in rcnormalon-subtracted 
schemes), the QCD potential calculated by lattice 
simulations, and the phenomenological potentials 
in this region. Combining these analyses, we 
find that all these potentials agree well with one 
another. It appears that the non-perturbative ef- 
fects, if they exist, are comparable in size to the 
perturbative uncertainty (apart from those in the 
r-independent part of the potentials, which has 
not been constrained by these analyses). 

Phenomenological potentials in the above 
range may be represented by a Coulomb-plus- 
lincar potential, which becomes steeper than the 
Coulomb potential at larger distances. We can 
understand why E tot (r) also becomes steeper 
than the Coulomb potential in perturbative QCD 
[4]. After realizing that the 0(A) renormalon 
should be cancelled, it is natural to define a strong 
coupling constant ap{p) from the interquark 
force: 

F(r)— IW)-- Cl,^. (2) 

Since the 0(A) renormalon in Vqcd(^) is r- 
independent, it is killed upon differentiation. 
ai?(l/r) grows at IR due to the running (the 
first two coefficients of the /3-function are uni- 
versal), which makes \F(r)\ stronger than the 
Coulomb force at large distances. This means 
that VqcdM, after subtraction of the renor- 
malon, becomes steeper than the Coulomb po- 
tential. 

2.2. Bottomonium spectrum 

The bottomonium energy levels can be com- 
puted in series expansions in as within the 
boundstate theory based on perturbative QCD. 
Presently the full corrections up to 0(a s rrib) are 
known [6]. It is illuminating to compare the se- 
ries expansions before and after the cancellation 
of the 0(A) renormalons. For instance: 
•T(1S) : 3 

M T(1S) = 9.94- 0.17- 0.20- 0.30 GeV (3) 
= 8.41 + 0.84 + 0.20 + 0.013 GeV (4) 

•T{2S) A 

M T(2S) = 9.94-0.10- 0.19 -0.45 GeV (5) 
= 8.41 + 1.46 + 0.093 + 0.009 GeV(6) 

3 The parameters are taken as fi = 2.49 GeV, as(/i) = 
0.274, m b = 4.20 GeV, mj, iPole = 4.97 GeV. 
4 The parameters are taken as ^ = 1.09 GeV, = 
0.433, m b = 4.20 GeV, m biPole = 4.97 GeV. 
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Fig. 2: Bottomonium spectrum up to n — 3 as 
calculated in various frameworks compared to the 
experimental data. 



The upper lines correspond to the pole-mass 
scheme and the lower lines to the MS-mass 
scheme. As can be seen, the convergence prop- 
erty improves dramatically when the MS mass is 
used instead of the pole mass. In this scheme, 
the perturbative series for the bottomonium en- 
ergy levels show converging behaviors up to the 
states with the principal quantum number n = 3 
[7,8]. 

Fig. 2 compares the bottomonium spectrum cal- 
culated in various frameworks with the exper- 
imental data: BSV(1161) and (1181) represent 
the fixed-order perturbative QCD predictions up 
to 0(a s m,b) for the input as(Mz) = 0.1161 and 
0.1181, respectively [8]; RS(1181) represents the 
perturbative QCD prediction which includes in 
addition part of the higher-order corrections, in 
particular the full (D(a s mb) corrections to the fine 
structure, for as(Mz) = 0.1181 and /i = 3 GeV 
[9] ; Eichten-Quigg represents the prediction of the 
phenomenological potential model [10]. The level 
of agreement between the perturbative prediction 
(RS) and the experimental data approximates to 
that of the Eichten-Quigg model. 

2.3. Physical picture [7] 

When the 0(A) renormalons are cancelled, the 
major part of the energy of a bottomonium state 
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can be written as 
E « 2m b + 



d 3 q 4-Ka s (q) 
(2^ ° F ^^^ 



(7) 



1 < q<m b 



where R represents the spacial size of the bound- 
state. It shows that the energy is mainly com- 
posed of (i) the MS masses of b and b, and (ii) 
the self-energies of b and 6 originating from glu- 
ons whose wavelengths are shorter than the size 
of the boundstate 1/rrib J$ A < R. As expected, 
contributions of IR gluons (A > R) have decou- 
pled from eq. (7) . The potential energy between b 
and b turns out to be much smaller than the self- 
energies (ii). This observation is reminiscent of 
the constituent quark mass picture for explain- 
ing the masses of light hadrons, if we regard 
the energy (ii) as the difference between (state- 
dependent) constituent quark masses and current 
quark masses. 

In Fig. 2 the level spacings between consecu- 
tive n's are almost constant, in contrast to the 
Coulomb spectrum in which the spacing decreases 
rapidly with n as 1/n 2 . This is because the sec- 
ond term of eq.(7) grows rapidly as the size of the 
boundstate increases. Namely, the rapid growth 
of as (q) at q ~ as R increases pushes up 
the energy levels of the excited states. This gives 
a microscopic description for how the interquark 
force between a color-singlet heavy quark pair be- 
comes strong at large distances (within the range 
where the perturbative prediction is still valid) . It 
is the rapid growth of the self-energies (ii) rather 
than of the potential energy that gives the domi- 
nant effect. 

3. Cancellation of C(A 3 ) Renormalons 

Let us briefly state a historical background re- 
lated to this subject. As already noted, the en- 
ergy of a static bb pair, eq.(l), contains 0(A 3 r 2 ) 
renormalon, where we may replace r by the typi- 
cal size R ~ {asmb)~ l of the boundstate. Within 
the potential-NRQCD framework, it was shown 
[11] that this renormalon can be absorbed into 
a non-local gluon condensate, i.e. the perturba- 
tive uncertainty can be factorized and replaced 
by a non-pcrturbative parameter of the same di- 
mension ~ A 3 . On the other hand, it has been 
known for a long time [12] that the leading non- 
perturbative corrections to the quarkonium en- 
ergy levels are 0(A 4 ) since they are proportional 
to the local gluon condensate (G^G^). Thus, 
there is an apparent mismatch in the power of A 



between the two quantities. 

It has recently been shown [13] that, if the off- 
shellness of b and b is incorporated properly, it 
provides an additional suppression factor of order 
A/(a s rrib) to the renormalons in the quarkonium 
energy levels: 



SE ~ A x 



A< 



{a s m b ) 2 



a s m,b 



(8) 



Hence, the dimension of the perturbative uncer- 
tainty becomes the same as that of the leading 
non-perturbative corrections (including the pow- 
ers of as and m& in the denominator). Moreover, 
convergence of the perturbative series improves 
by this effect. 

Intuitively the suppression mechanism may be 
understood as follows. Large offshellness of b and 
b corresponds to short rescattering time of b and 
b inside the boundstate. If the rescattering time 
At ~ (a|m b ) _1 is shorter than the hadronization 
time ~ A -1 , b and b will get distorted before IR 
gluons surround them and an energy is accumu- 
lated by antiscrecning effects. Thus, we expect 
the offshellness to act as an infrared cutoff to the 
effects induced by gluons' timelike propagation. 

In the (leading) kinematical configuration rel- 
evant to formation of the bottomonium states, 
gluons exchanged between b and b have momenta 
q^ where <C \q\. Therefore, in the conven- 
tional approach we take the instantaneous gluon- 
exchange as the leading-order and incorporate 
perturbations to it. On the other hand, a dif- 
ferent kinematical region, — q 2 ~ A 2 , contributes 
to renormalons. Because of this, in order to esti- 
mate renormalons in the bottomonium spectrum 
accurately, we need to make an approximation 
valid in both of the above kinematical regions. 

An analysis following this observation reveals 
that there is a rich mathematical structure un- 
derlying the renormalon suppression effect. For 
instance, the expansion about the instantaneous 
potential changes the analyticity of the Borcl 
transform E[u] of the bottomonium energy level. 
IR renormalons, defined as poles of E[u] in the it- 
plane, arc located at u = 2,3,4,.... But once 
the expansion is made, poles are generated at 
u=|,|,|,..., while at the leading-order of the 
expansion, the original poles at u = 2,3,4,... 
disappear. For example, at it = |, there is a can- 
cellation of poles in the combination 



E\u] 



A 3 " 



2u 



it -3/2 ' 



where A = \C F a S - (9) 
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A represents the offshellness or non-instantaneity 
of the potential. However, if A is considered to 
be small and an expansion A 3 ~ 2u = A 3 x (1 — 
2u log A + . . .) is performed, there remains an un- 
cancelled pole at u = §. Hence, the perturbative 
uncertainty becomes larger by this expansion. 

4. Conclusions 

There has been important progress in the the- 
ory of heavy quarkonium states, which extended 
the predictive power of perturbative QCD beyond 
what could be achieved before. Perturbative pre- 
dictions can be made accurate once IR degrees of 
freedom are eliminated properly from the calcu- 
lations. Renormalon dominance picture suggests 
that uncertainties of the perturbative prediction 
for the bottomonium spectrum are of moderate 
size, of C(A 3 /(a|m^)) if the short-distance mass 
is used, of C(A 4 /(a|m 3 )) if the offshell effects are 
further incorporated. Evidence has been found 
for the static QCD potential and the bottomo- 
nium spectrum, which supports the hypothesis 
that the perturbative predictions agree with the 
full QCD predictions within the estimated mod- 
erate uncertainties. 

Based on the perturbative QCD predictions, we 
obtained new microscopic pictures on the compo- 
sition of the heavy quarkonium spectrum and on 
the behavior of the interquark force in the inter- 
mediate distance region. Furthermore, the quark 
offshell effect as an IR cutoff in timelike processes 
was recognized, whose nature is essentially non- 
local in time. 

So far, most important applications of the 
progress are the precise determinations of the MS 
masses of b and c quarks: e.g. nib = 4190 ± 
32 MeV, m c = 1237 ± 60 MeV have been ob- 
tained based on the above hypothesis that non- 
pcrturbative corrections can be absorbed into 
perturbative uncertainties [8]. 
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